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We have developed a scheme of Gor’kov Green’s functions to treat the impurity effects on the 
Caroli-de Gennes-Matricon (CdGM) mode in superconductors (SCs) by improving the Kopnin- 
Kravtsov scheme with respect to the coherence factors and applicability to various SCs. We can 
study the impurity effects while keeping the discreteness of the energy spectrum, in contrast to the 
quasiclassical theory. We can thus apply this scheme to SCs with a small quasiclassical parameter 
fcp£o (which is the product of the Fermi wavenumber kp and the coherence length £o in a pure SC at 
zero temperature) and/or in the superclean regime A m i n ir 1 (A m i n i and r denote, respectively, 
the level spacing of the CdGM mode called the minigap and the relaxation time for the CdGM mode 
and we take h = 1). We investigate the impurity effects as a white noise for a vortex in an s-wave 
SC and two types of vortex in a chiral p-wave SC for various values of the quasiclassical parameter 
and impurity strength (from the moderately clean regime to the superclean regime) and confirm the 
validity of this scheme. 


I. INTRODUCTION 

In a vortex core in a superconductor (SC), there are 
electronic bound states called the Caroli-de Gennes- 
Matricon (CdGM) mode, which has been predicted 
by Caroli et alX and observed experimentally^ and 
numerically^. These excitations are described as the so¬ 
lutions to the Bogoliubov-de Gennes (BdG) equation. 
The CdGM mode makes an important contribution to 
phenomena in low-energy physics in vortex states, such 
as flux flow states^. The flux flow conductivity in a 
dirty SC is often explained by Bardeen-Stephen theory^, 
in which the vortex core region is treated as a normal 
state. In a clean SC, however, the impurity scattering of 
the CdGM mode is the main mechanism of dissipation in 
the flux flow states. Moreover, particle-hole asymmetry 
in the CdGM mode causes a flux flow Hall effect. The 
impurity effects on the CdGM mode play a key role in 
transport properties in mixed states. 

Recently this vortex core mode has attracted much at¬ 
tention owing to the existence of Majorana quasiparticles 
in zero-energy states in topological SCs££. The Majo¬ 
rana quasiparticle is its own antiparticle and is expected 
to exist at the edge or at topological defects such as a vor¬ 
tex core of a topological SC. An idea for the realization of 
topological quantum computation (TQC)^ has been pro¬ 
posed, which makes use of non-Abelian statistics in the 
Majorana zero-energy states in vortex cores. However, it 
requires adiabatic braiding operations of Majorana zero- 
energy states. This adiabatic time scale should be less 
than the inverse of the minigap A m ; n ;, which is the level 
spacing of the CdGM mode in the vicinity of zero energy 
(in this paper we take h = 1). Typically the minigap is 
on the order of A m ; n i ~ Ao/(fcp£o) where Ao, kp, and 
£o are, respectively, the modulus of the pair potential in 
the bulk far from the vortex core, the Fermi wavenum¬ 


ber, and the coherence length at zero temperature. In 
Sr 2 Ru 04 , which is expected to be a chiral p-wave SC, 
A m ini is estimated to be on the order of lmK—. The 
magnitude of the minigap is relevant for the realization 
of physics associated with Majorana states. Since real de¬ 
vices include impurities, it is an important issue how ran¬ 
domness affects the minigap and the possibility of TQC. 

There are three regimes for the inverse of the relax¬ 
ation time t“. (i) the dirty regime t” 1 > Abuik, (ii) the 
moderately clean regime A mini < r” 1 <C Ab u ik, (hi) the 
superclean regime rff 1 -C A m mj. It is reported that im¬ 
purity effects depend on the number of impurities -/V; mp 
in a vortex core as well as the relaxation time r n , partic¬ 
ularly in the superclean regime. In other words, r n can 
be expressed as A m i n ir n oc (-/Vj mp 0 2 ) -1 , where 6 is the 
Born parameter and characterizes the type of scatterers, 
and the magnitudes of both N[ mp and 6 are crucial for 
determining impurity effects, especially for large r n . 

In the presence of a few impurities with the Born pa¬ 
rameter being not small, 1 /v^f? <C 9 < 1, the 2A m ; n i- 
periodic structure of the energy spectrum appears rather 
than the A m ; n i-periodic oneiir— . These kinds of impu¬ 
rity always give rise to level mixing, and the Landau- 
Zener transition causes the dominant absorption in flux 
flow states in the (moderately clean and) superclean 
regimes. For an unconventional superconductor such as 
a chiral p-wave SC or d-wave SC, these impurities lead 
to the density of states (DOS) around the Fermi energy, 
corresponding to the midgap bound states^. 

The adiabatic dynamics in TQC requires, however, ex¬ 
ceptionally small absorption, making it also important to 
consider the nondissipative regime for vanishingly small 
0. We assume there are weak but many scatterers in 
a vortex core. This kind of impurity is called white 
noise disorder. In the dirty regime and moderately clean 
regime in which the Born parameter is very small, analy- 
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sis with random matrix theory (RMT) has clarified that 
the level statistics obey class C of the Altland-Zirnbauer 
classification— -il. Impurity scattering rates have been 
obtained by analysis using the quasiclassical theory!**— 
In these studies, spectra form peaks with widths. The 
analysis by the quasiclassical theory has also clarified the 
difference in the types of SC, such as s-wave and chiral p- 
wave SCs. The difference stems from coherence factors, 
and the cancellation of coherence factors corresponding 
to Anderson’s theorem in the vortex core occurs in one 
of the vortices in a chiral p-wave SC— 

For the adiabaticity in TQC, we consider the super¬ 
clean regime with white noise disorder. In this regime, 
we need to treat discrete levels, hence the quasiclassi¬ 
cal theory is not applicable. Kopnin and Kravtsov have 
developed a scheme of Gor’kov Green’s functions with 
an impurity self-energy to investigate the flux flow Hall 
effect in a three-dimensional type-II s-wave SC2£, where 
they have treated the vortex core states while keeping the 
discreteness of their levels. The form of the coherence fac¬ 
tors and thus their way of application to other SCs are, 
however, not clear in their scheme. In this paper, to over¬ 
come these difficulties, we derive a suitable scheme for 
Gor’kov Green’s functions with an impurity self-energy 
in two dimensions and apply it to a chiral p-wave SC 
with two types of vortex as well as an s-wave SC with a 
single vortex. This scheme is easily extended to three- 
dimensional cases if we do not consider numerical costs. 
We confirm the validity of this scheme by comparing the 
results for impurity effects in a system with o = 100 
with those of the quasiclassical theory. We also investi¬ 
gate the impurity effects in the large-A m j n iT regime (r 
denotes the relaxation time for the CdGM mode), where 
the system is not in the quasiclassical regime, and find 
that the impurities affect spectra in the low-energy re¬ 
gion in a way similar to those in the quasiclassical limit, 
except for the zero-energy state. 

In the following, we use the unit /cb = 1 in addition to 
H = 1 for convenience. 


II. METHOD 


energy. We define Green’s functions by 


G(xi,x 2 ) 


( G(x i,x 2 ) F(x 1 ,x 2 ) 

y-F^(x 1 ,x 2 ) G(xi, x 2 ) 


( 2 ) 


Each element in Eq. m represents a 2x2 matrix, the 
matrix elements of which are given by 

(G(xi,x 2 )) a p = (Tijj a (xi)ilj^(x 2 )), (3) 

(G{x 1 ,x 2 )) a p = —(Tipl e (xi)'ipp(x 2 )), (4) 

(F( Xl,X 2 )) a p = (Tip a (xi)ipp(x 2 )}, (5) 

{F\xi, X 2 )) a /3 = (T'lf)l l (xi)lpp(x 2 )). ( 6 ) 

The symbols a and /? denote the electron spins f or 4~ 
We introduce the Fourier transformation for relative time 
f = n - t 2 : 


G(x i,x 2 ) = G(ri,r 2 ;r) = T^G(r i,r 2 ;v n )e lu>nT , 

n 

( 7 ) 

where uj n denotes the fermionic Matsubara frequency. 
Within the self-consistent Born approximation (SCBA), 
Ei m p is described as 

F n „ 

ShnpO; w n ) = —^G(r, r; u n ) = 7 G(r, r; w„). (8) 

ttAG 

N n denotes the density of states per spin at the Fermi 
energy ep in the normal state. T n is the impurity scatter¬ 
ing rate in the normal state and related to the relaxation 
time T n as r n = l/r n . G^ can be obtained from the 
BdG equation without impurities: 

^BdG( r ) M ^( r ) = E K u K (r). (9) 

E^BdG ( r ) th e BdG Hamiltonian for the 

pure SC. The basis takes the form Ux{r) = 
t {u K ^{r),UKi{r) 1 v K ^{r),VKi{r)) with index K. The 
expression for G^ is given by 


A. General approach 

We consider the Gor’kov equation for superconductors 
with a single vortex and rotational symmetry. The sys¬ 
tem has a disc geometry of radius R from the vortex 
core at the origin. The Gor’kov equation in real space is 
transformed into the Dyson equation in real space: 


G (0) (r , r'\ u> n ) 


K 


z u K {r)(u K (r')y 

E k - i iO n 


( 10 ) 


Since the system without impurities has rotational sym¬ 
metry, the modified angular momentum J z , 


Jz — L z 1 + lia z f z + 1 2 t z , 


( 11 ) 


G(r,r';w„) = G (0) (r, r'\u n ) 

+ J dn& 0) (r i ; tdn)Dj m p(7*1, Cd n )G(Vi , V ,ld n ), 

(1) 

where " denotes a 4 x 4 matrix, G(G^) Green’s function 
with (without) impurities, and Ei mp an impurity self- 


commutes with the Hamiltonian and 

its eigenvalue l is a good quantum num¬ 
ber. We define a unitary operator Ui{&) = 

diag(e i( ' _Zl_i2)e , e i(i+Zl_i2)e , e i( - l+h+l2 ) B , e i ( l - ll+l2 ' >B )/V2ir, 
and describe the eigenvector Uk{t) as 

u K (r) = Ui(9)n^(r). (12) 
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Here we introduce (r, 8) as two-dimensional polar coordi¬ 
nates. The basis is described as (e r ,eg) = R(9)(e x ,e y ), 
where Rid) is a 2 x 2 matrix that rotates the vector 
around the z-axis by angle 6. The index K is rewritten 
as the set of the angular momentum l and the quantum 
number of the radial part u. To eliminate the 9 depen¬ 
dence from the BdG equation, we multiply both sides of 
Eq. © by U}'(9) from the left and integrate them with 
respect to 9 over 0 to 2i r. The BdG equation for this 
radial part 7 Z-i tI/ (r) is obtained as 

uf\r)n liV {r)=E l ^ v {r), (13) 


where cldU^ (9)H^J G (r)Ui(9) = 5i^Ti[ 0 \r). Green’s 
function can be written as 


G (0) (r,r'; 
G, (0) (r, r'\ 


l 

, Y- f z 'R.i, v (r){Ki, v (r')y 

Eiv — i io n 

IS 


(14) 

(15) 


Substituting Eqs. (fl8l) and (1201) into Eq. we obtain 
an equation for a: 


&l,vi ,V 2 (i^n) — 2^ 


E 

' y. y„ 


M, 


w 


W [(l-^'(i Wn))~ 


(iw„ - E Vv d)(vjj n - Ei V2 ) 


( 21 ) 


where 
M 11 ' , , 

Vl V <2 V ' 2 


r <M' R -/U (r)] [7L\, . {r)f z n l}U2 (r)], 

( 22 ) 


and these matrix elements [Eq. (l22l) ] represent coherence 
factors. The density of states can be evaluated from the 
Green’s functions obtained above as follows: 


N(u>) = J drlm 


- Trf z G(r,r;w n ) |. 




= E Im 

l,IS 


- [(1 - ai(u + i5)) x ] 
7r(w + iS - Ei v ) 


E^M- 

l,is 

(23) 


Although impurities break the rotational symmetry, scat¬ 
tering processes do not change the angular momentum af¬ 
ter averaging over the configurations of impurities. This 
means that Green’s function G also takes the following 
form: 

G(r,r';u, n ) = J^I7 l (9)G l (r,r';u, n )tf l f (9'), (16) 

l 

and we can obtain the radial part of the Dyson equation 
as 

G;(r, r';u„) = G[ 0) (r, r'\ u n ) 

+ / r 1 dr 1 Gf\r,r 1 \uj n )% m - p (ri]UJ n )Gi(r l ,r'-,u> n ), 

(17) 

^imp (^li ^n) — / J Gl(v\, T\ : iO n ). ( 18 ) 

l 


B. Single-mode approximation 

We are interested in the impurity effects on the CdGM 
mode and assume that the main contribution to this 
mode is from the CdGM mode itself. We thus restrict 
v to that for the CdGM mode, to which we assign 
v = 0. We will examine the validity of this restriction in 
Sect. IIII Cl We redefine a and M as follows: 

cn(iu n ) = (E t ,o - iu n )ai$fi(\u n ), (24) 

Mi,v = My o, 0 , 0 - (25) 

Equations m and (1^1) are then reduced to 

= f E -W - M r\ ■ • ( 26 > 

27 T E vo - cri>(vjj n ) - 1 UJ n 

My = J rdr\n[ 0 (r)f z n^ 0 (r)\ 2 , (27) 


We define the scalar quantities 


&l,is,is' (l^n) — 


rdrK] M (r)- 


J imp 


(r-,u> n )f z 1Ziy ( r) 


Eli/' kw 


(19) 


and they constitute a matrix di(iw n ): (d’i(iu> n )) L ,y = 
ai, v y( iw n ). With this modified self-energy tf, the Dyson 
equation can be solved formally, and Green’s function is 
obtained as 


Gi(r, r';u)n) = ^ 

isis' 


f z T^Uis{r) 
Ehs 1 


E (< J l(. iuJ n )) k 


K-l ,(/) 


isis' 


, Eh, — i u n 


di(iuj n )) % ] vv ,' R -\ y (r'). 


( 20 ) 


and Green’s function Gi is rewritten as 


G z (r, r';uj n ) 


f z n l}0 (r)(n L0 (r')Y 
Ei ,o - cq(iw n ) - i uj n ' 


(28) 


Hereafter, we call this scheme the improved Kopnin- 
Kravtsov (iKK) scheme. The iKK scheme is advanta¬ 
geous to the scheme derived by Kopnin and Kravtsov^; 
Eq. (l26l) has the coherence factor Miy, which was not 
clear in Eq. (1) of Ref. [2C]. The coherence factor de¬ 
pends on the pairing symmetry, and the iKK scheme can 
thus reflect differences among various types of SC. Since 
Eqs. (l26l) (|28l) are constructed with the eigensystems of 
the BdG equation formally without reference to their an¬ 
alytic expressions, we can easily apply this method to 
other superconducting systems with a single vortex for 
which an analytical solution to the CdGM mode is not 
available. 
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FIG. 1. (Color online) Impurity scattering rate F as a func¬ 
tion of energy eigenvalue E for L z = 0, 1, 2. QC means the 
quasiclassical calculation and solid lines are calculated by the 
quasiclassical theory^. Except for the case with L z = 0, the 
results of the iKK scheme are in good agreement with those 
of the quasiclassical theory. 


III. NUMERICAL STUDY 

We consider two cases, where the quasiclassical pa¬ 
rameter /cf£o is 100 and 5. We calculate the DOS within 
the single-mode approximation introduced in Sect. Ill Bl 
We take the units of energy and length for Ao and £o, 
which are the modulus of the pair potential in the bulk 
far from the vortex core and the coherence length defined 
as £o = tr/Ao, respectively. In this paper, we consider a 
single vortex in an s-wave SC and two types of vortex of 
chiral p-wave SCs. In chiral p-wave SCs, there are two 
types of vortex, each of which cannot be transformed into 
the other by any symmetry operation. This kind of in¬ 
equivalence is inherent to SCs in which the time-reversal 
symmetry is broken even in the absence of a magnetic 
field; recall that a vortex and antivortex in an s-wave SC 
are transformed to each other via time-reversal operation. 
Each type of vortex in chiral p-wave SCs is specified by 
the relative sign of the vorticity and chirality. 

For s-wave and chiral p-wave SCs, the 4x4 BdG Hamil¬ 
tonians can be reduced to the spinless representation. 
The explicit forms are given by 


We assume the spatial profile of the pair potential to be 

A(r) = A 0 tanh(r/£ 0 )e ±ie . (31) 

We take the sign of the chirality in the chiral p-wave 
SC in Eq. (l30l) as positive. For the chiral p-wave SC, 
the positive sign in Eq. m means that the vorticity is 
parallel to the chirality and the negative sign means that 
the vorticity antiparallcl to the chirality. The former is 
labeled as L z = 2 and the latter as L z = 0, and we call the 
vortex with L z = 2 ( L z = 0) the “parallel (antiparallel) 
vortex”. 

In order to obtain eigenvalues and eigenfunctions, we 
use the Fourier-Bessel expansio n 2 — 22 . We use the radius 
R = 10£o and two cutoffs, l c = 100, 30, and N c = 200 
for l and is, respectively. N c denotes the number of zero 
points of the Bessel functions we have used. In the follow¬ 
ing calculations, we consider T = 0 and replace iui n with 
w + i(5 in the above scheme represented by Eqs. (EHD (EHD- 
8 is infinitesimal and set to 10 _5 A o . Impurity effects are 
evaluated through the width of the spectrum, denoted 
as F. We define T as the half width at half maximum, 
which we also call the impurity scattering rate. The re¬ 
laxation time r for the CdGM mode is associated with T 
as r ~ i/r. 

We remark on how to construct the Majorana state 
localized around the vortex core in numerical calcula¬ 
tions for finite-size systems. In chiral p-wave SCs with 
an infinite system size, there are two zero-energy states 
with l = 0. Their energies become finite but sufficiently 
small and opposite in sign to each other when the sys¬ 
tem size is finite and sufficiently large. We use the la¬ 
bels (l = 0 , is = +) for the positive-energy state and 
(l = 0, is = —) for the negative-energy state. Moreover, 
both their wave functions simultaneously have weight in 
the vortex core and near the edge. In order to treat these 
states within the single-mode approximation, we define 
the energy and the state labeled by is = 0 and l = 0 
as £ 0 ,o = (E 0j+ ± E 0 -)/2 and = [R- Q ,+ (r) ± 

Ko,-(r)\/V2 t respectively. We choose the sign of ± so 
as to obtain the bound state in the vortex core. 


A. Results: = 100 


/-V 2 

—- n D(r) 

DM 5 -+/ 


u K {r) 

v K {r) 


= E k 


u K (r) 

v K (r) 


(29) 


when 


.21 


s-wave, 


D(r) = < 


'A(r) = (D(r))* 
t[ 5±A(r)] + A(r)S±| = - (-D(r))* 
_ d± = d x ± i d y chiral p-wave. 


(30) 


For this quasiclassical parameter, we expect that this 
method will reproduce the result calculated by the quasi¬ 
classical theory^. The normal scattering rate T n can be 
changed as a parameter. First, we set T n = 10 _2 7 tA o for 
an s-wave SC ( L z = 1) and a chiral p-wave SC with 
L z = 2, and r n = 10 _1 7 tAo for a chiral p-wave SC 
with L z = 0. In these cases, the widths T are summa¬ 
rized in Fig. Q] which are evaluated from the spectra in 
Figs.EKa)-(f). Figures[2ja)-(c) show the DOS calculated 
in the iKK scheme for each angular momentum l , which 
takes discrete values. Figures 2(d)-(f) show the angular- 
resolved local density of states (LDOS)^ calculated by 
the quasiclassical theory for some impact parameters b 
(b/t; o = O.Oln, n = 0,1,2, ■ ■ ■), which can take arbitrary 
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FIG. 2. DOS for CdGM modes of the s-wave SC ((a) (d), L z = 1) and the chiral p-wave SCs with the parallel vortex ((b) (e), 



FIG. 3. (Color online) Impurity effects for the normal scattering rate r n = 10” 4 7rAo and the quasiclassical parameter /cf£o = 
100. (a) and (b) show the DOS for the s-wave SC and the chiral p-wave SC with the parallel vortex, respectively. The widths 
of the spectra for L z — 0,1, 2 are shown in (c). 


values. In Figs. E](d)-(f), we have introduced another 
coordinate system spanned by (e s ,e{,) = R(a)(e Xl e y ), 
where a is the angle between kp and e x . We assume ro¬ 
tational symmetry; thus, the angular-resolved LDOS is 
independent of the value of a for the fixed set of (s, b). 
We calculate the LDOS at s = 0, referring to Fig. 3 
in Ref. QiJ We remark on the relationship between the 
DOS and LDOS. In these cases, the angular momentum 
is approximately described by the impact parameter as 
l = —kpb, which is not quantized. Because of this re¬ 
lationship, the scattering rate and peak energy for the 
state labeled by l correspond to those for the state la¬ 
beled by (kp,b). Therefore, the DOS and LDOS have 
similar structures since the scattering rate and peak en¬ 


ergy mainly characterize the shape of the spectrum. 

Figure 0(a) shows the energy dependence of the DOS 
Ni t0 (E) [defined as Eq. (l23l) ] for the s-wave SC. This spec¬ 
trum is very similar to that obtained by the quasiclassical 
theory (Fig. Od))- We find that the spectra also overlap 
their neighbor states regardless of their discreteness. In 
this parameter region, we can thus consider the system to 
be in the quasiclassical regime. In fact, the typical width 
of the spectrum T ~ T n = 10~ 2 7rAo is compatible with 
the level spacing, A 0 /(/cf^o) ~ 10“ 2 Ao; thus, the levels 
are regarded as continuous. The DOS and LDOS of the 
chiral p-wave SC with the parallel vortex (L z = 2) are 
shown in Figs. Ob) andO e )- hr this case, each spectrum 
overlaps its neighbors enough for the system to be in the 
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FIG. 4. (Color online) Impurity effects for the normal scattering rate T n = 10 _2 7rAo and the quasiclassical parameter fcF^o = 5. 
(a) and (b) show the DOS for the s-wave SC and the chiral p-wave SC with the parallel vortex, respectively. The widths of the 
spectra for L z = 0, 1, 2 are shown in (c). 


quasiclassical regime, and is in good agreement with that 
obtained by the quasiclassical theory as for the case of 
the s-wave SC. The energy dependence of the spectrum 
broadening is similar to that for the s-wave SC, where the 
spectrum with lower energy in the CdGM mode broadens 
more. 

For the chiral p-wave SC with the antiparallel vortex 
(L z = 0), the spectra are sharp, as shown in Fig. (Sic), 
i.e., the spectra are hardly affected by impurities. These 
small scattering rates are qualitatively the same as those 
obtained by the quasiclassical theory. Because of this ro¬ 
bustness, there are almost no overlaps among the spectra 
shown in Fig. [2j)c) , while in the quasiclassical theory the 
small overlaps stem from the continuous spectra shown 
in Fig. [2]T). The widths obtained by this method are 
slightly different from those obtained by the quasiclassi¬ 
cal theory, as can be seen from the blue solid line and blue 
triangles in Fig. [Q For these values of o and r n /A 0 , 
the system with L z = 0 can no longer be considered to 
be in the quasiclassical regime. This implies that we can 
consider the system is in the quasiclassical regime when 
there are sufficiently large overlaps among the spectra. 
In contrast to the s-wave SC and the chiral p-wave SC 
with L z = 2 mentioned above, in terms of the energy 
dependence of scattering rates, spectra around zero en¬ 
ergy are more robust against impurities than those at 
higher energies in CdGM mode. This behavior is quali¬ 
tatively the same as the quasiclassical behavior. We thus 
see that the energy dependences of the scattering rates 
of CdGM modes are similar for the quasiclassical and 
quantum (A m i n ;r ^1) limits. 

From these figures (Figs. Ufa)-(f)), we thus find that 
the result obtained by the iKK scheme and the quasiclas¬ 
sical theory are in good agreement except for the chiral 
p-wave SC with L z = 0; Figures He) and Hf) exhibit 
a small discrepancy with respect to the overlaps of the 
spectra. The system with L z = 0 is not in the quasi¬ 
classical regime even for fcp£o = 100, as we mentioned 
above. 

Next, we consider another case, in which the systems 
are very clean, and hence the widths are small and the 
overlaps among the spectra are also small. We consider 


only the s-wave SC and the chiral p-wave SC with L z = 2 
and set T n = 10 _4 7 tAo. Both of these SCs do not have 
overlaps of the spectra (as shown in Figs.H a ) andHb)). 
Actually, the widths T of the s-wave SC and the chi¬ 
ral p-wave SC with L z = 2 are on the order of 10” 3 Ao 
(Fig. H c ))> while the level spacings are on the order of 
10 _2 Ao- The tendency of the impurity effects is the same 
as in the above results and quasiclassical results. The 
chiral p-wave SC with L z = 2 has larger impurity scat¬ 
tering rates at lower energies but the zero-energy state 
has an exceptionally small scattering rate. This is be¬ 
cause the matrix element M 0 ,o [Eq. (1271) ] is zero due to 
the character of the zero-energy state. When the spectra 
are discrete and the contributions from neighbor states 
are relatively small, this character appears, in contrast 
to the above state with L z = 2. In this sense, the vortex 
type of the chiral p-wave SC with L z = 2 is different from 
that of the s-wave SC. 


B. Results: /cf£o = 5 

The system with fcp£o = 5 is not in the quasiclassical 
regime because of the large level spacing. Let us discuss 
some differences between the results for the small qua¬ 
siclassical parameter (/cf£o = 5) and large quasiclassical 
parameter (&f£o = 100). Although the normal scatter¬ 
ing rates (r n = 10~ 2 7tAo) are the same in these two 
cases, the resultant scattering rates V in Fig. H c ) are 
larger than those in Fig. |T] We infer that a large level 
spacing leads to a smaller level repulsion, which yields 
the wider spectra. The impurity effects thus depend on 
/cf£o- This dependence is not accessible in the quasi¬ 
classical theory because the conventional quasiclassical 
theory neglects the contributions from O(l/(fcp£o))- In 
Fig. He), the scattering rate T is small only for the zero- 
energy state in the chiral p-wave SC with L z = 2, as 
well as that in Fig. He). Overall behaviors are qualita¬ 
tively the same as those in Fig. H c )- We note that since 
kpVpT = /cpCo^o/Fn = 500/7T 1, weak-localization 

effects are small, regardless of the small quasiclassical 
parameter. 
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FIG. 5. (Color online) ^-dependences of the widths of the 
states (a) l = 0 and (b) l = 1 in the CdGM mode for the cutoff 
l c = 3. N is the number of modes, counting from the lowest- 
energy mode, in the extended states with positive energies. 
The type of the vortex is L z = 2. Solid squares denote the 
width with N modes. The dashed lines and solid lines denote 
the width within the single-mode approximation for l c = 3 
and lc = 30, respectively, for reference. 
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FIG. 6. (Color online) ^-dependences of the widths of the 
states (a) l = 0, (b) l = 1, and (c) l = 2 in the CdGM mode 
for the cutoff l c = 5. The symbols are the same as in Fig. [5] 
except that the dashed lines denote the width for l c = 5. 


C. Contributions from the modes in the bulk 

In this section, we evaluate the contributions from the 
extended states in the bulk with energies above Ao using 
Eqs. m ©• We consider the chiral p-wave SC with 
L. = 2 and use o = 5 and T n = 10 _2 7rA o . To reduce 
numerical costs, we set the cutoffs l c to 3 and 5, where l c 
gives the maximum of |Z|. For a given l, we take account 
of 2 TV + 1 modes, which consist of TV modes with posi¬ 
tive energies, N modes with negative energies, and one 


CdGM mode. The results for l c = 3 (Fig. 0 are almost 
the same as those for l c = 5 (Fig. 0. In these figures, 
open squares represent widths T j(JV) (i = 0,1,2: the la¬ 
bels of excited states in the vortex core) and blue solid 
lines represent the widths for TV = 0 and l c = 30. Here, T 
used in the above sections corresponds to r^JV = 0). As 
TV increases, the widths I 1 * (TV) appear to converge. The 
contributions to the width of the zero-energy state are 
about 10% of r 0 (TV = 0) in Figs. 0 a) and 0 a). More¬ 
over, the contributions to the widths of the first and sec¬ 
ond excited states with finite excitation energies are less 
than a few percents of Ti(TV = 0) and T 2 (TV = 0), respec¬ 
tively, in Figs. 0b), [6(b), and [He). There is little effect 
on excited states; thus, the single-mode approximations 
are sufficiently valid to obtain the scattering rates. 


IV. SUMMARY AND DISCUSSION 

We have developed a scheme based on the Gor’kov 
Green’s function in order to treat the impurity effects on 
bound states within the self-consistent Born approxima¬ 
tion (iKK scheme). This scheme properly takes account 
of the coherence factor and thereby it is applicable to 
various SCs, in contrast to the Kopnin-Kravtsov scheme. 
We have confirmed the validity of the iKK scheme by re¬ 
producing the result obtained by the quasiclassical the¬ 
ory for an s-wave SC and a chiral p-wave SC with L z = 2 
when the system is in the quasiclassical regime. For a 
chiral p-wave SC with L z = 0, on the other hand, the 
scattering rates calculated in the scheme are small and 
qualitatively the same but not quantitatively the same 
as those obtained by the quasiclassical theory because 
the spectra have almost no overlaps and are regarded as 
discrete; T ~ 1/r < A m ; ni . 

This scheme has three advantages over the quasiclas¬ 
sical theory. First, the scheme is applicable to the cases 
where o is small and the quasiclassical theory is in¬ 
valid. In this parameter region, the impurity scatter¬ 
ing rates calculated by the scheme are quantitatively dif¬ 
ferent from those calculated by the quasiclassical the¬ 
ory. Considering the difference in the case of small /cf£o 
and the small discrepancy for L z = 0 mentioned above, 
we conclude that the quasiclassical theory is valid when 
A mini r < 1. 

Second, in contrast to the quasiclassical theory, the 
iKK scheme is effective for calculating the impurity ef¬ 
fects on the zero-energy state of a topological SC such as 
the chiral p-wave SC in this paper. As a result, for a par¬ 
allel vortex in a chiral p-wave SC, the width of the zero- 
energy state is finite but exceptionally small compared 
with the widths of other excited states when we cannot 
regard spectra as continuous. (This property was not 
known in earlier analysis using quasiclassical theory^.) 
The width of the first excited state rather than the small 
width of the Majorana state is thus important in evalu¬ 
ating the minigap. Afo,o = 0 leads to this exceptionally 
small width when the spectra are discrete, as mentioned 


















in Sect. IIII Al Wc infer that this is associated with the 
character of the Majorana state. Although we can dis¬ 
cuss the impurity effects on the minigap and adiabaticity 
of Majorana states in the vortex core, there remains the 
issue of the origin of the small width of the zero-energy 
state. Even though we consider the extended states in the 
bulk, the width of the zero-energy state remains finite 
within the self-consistent Born approximation for ran¬ 
dom averaged impurities. It is, thus, not the single-mode 
approximation that causes the finite width of the zero- 
energy state. 

Third, the scheme will be useful for calculating the im¬ 
purity effects on the vortex core states in undoped Dirac 
SCs or SCs for which the derivation of the quasiclassi- 
cal equation is involved owing to the strong spin orbit 
coupling or Zeeman energy, provided the system has ro¬ 
tational symmetry. 

Lastly, we briefly remark on a possible future work on 
the two types of vortex in a chiral p-wave SC in terms of 


topological classification^!. A chiral p-wave SC belongs 
to the BdG class with broken time-reversal symmetry, the 
class D, which is one of the ten symmetry classes given 
by Altland and Zirnbauer, obtained by using the random 
matrix classification!!. We find, however, that the two 
types of single vortex in a chiral p-wave SC exhibit dif¬ 
ferent impurity effects. Our findings strongly suggest an¬ 
other possible classification in terms of symmetry and/or 
topology in the real space. 
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